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Abstract 

We investigate tlie electron addition spectrum in a class of Hubbard-like 

models which describe arrays of coupled quantum dots. Interdot tunneling 

leads to a sequence of two phase transitions separating a region of collective 

Coulomb blockade from a region where the Coulomb blockade of individual 

dots is maintained and a region where the Coulomb blockade is destroyed 

altogether. Observable experimental consequences of our theory are discussed. 
PACS numbers: 73.20.Dx, 73.40.Gk, 71.27.+a 



1 



Typeset using REVTgX 



Arrays of coupled quantum dots provide a novel system in which one can study 
"solid-state" physics on a much lower energy scale. Because of the dominant role played 
by electron-electron interactions in quantum dots, as manifested in the phenomenon of 
Coulomb blockade 0, it is expected that quantum dot arrays should mimic the physics of 
strongly-correlated lattice systems, and should therefore provide rather good experimental 
testing grounds for studying the predictions of Hubbard-type strongly-correlated electron 
interaction models. To date, theoretical treatments of quantum dot arrays have focused 
largely on arrays of metallic dots [Q, where the discreteness of the dot energy levels can 
be neglected, or arrays with nearly transparent barriers between dots 0|, where Coulomb 
blockade is absent. In this Letter, we investigate the electron addition spectrum in arrays 
of coupled semiconductor quantum dots in the Coulomb blockade regime and its vicinity. 
Our results directly reflect the intricacies of Mott-Hubbard correlations, and, in particular, 
show that various interaction-driven phase transitions should be observable in suitably 
fabricated quantum dot arrays. 

Motivated by the success in describing Coulomb blockade in a single quantum dot using 



an Anderson- type Hamiltonian |T^, we employ a generalized Hubbard model to account for 
the effects of quantum confinement, intradot Coulomb interactions, and interdot tunneling 
in a quantum dot array. We find a remarkably rich phase diagram for quantum dot arrays 
which generically exhibits three phases: for weak interdot tunneling, the quantum dot states 
are split into minibands but the Coulomb blockade of individual dots is maintained; for 
intermediate tunneling, the Coulomb blockade of individual dots is destroyed but there 
remains an energy gap due to interdot correlations (collective Coulomb blockade); and for 
strong tunneling, the Coulomb blockade is destroyed altogether. Our results provide new 
insight into an early experiment on transport through a one-dimensional (ID) array of 
quantum dots 0, as well as a possible explanation for anamolous degeneracies observed in 
the addition spectrum of a disordered quantum dot with several local minima . 

The system we wish to model consists of a linear or square array of identical quantum 
dots of the type described in Ref. [l^, separated from a metallic backgate by a thin insulating 



layer. We represent the Coulomb interaction between electrons on the same dot by an energy 
U = e^/C, the capacitive charging energy of the dot 0, and assume that interdot Coulomb 
interactions are screened by the backgate. The single-particle energy levels in the confining 
potential of an isolated dot are denoted by e^, a = 1, . . . , M (a specifies the quantum state 
of both the orbital and spin degrees of freedom), which we take to be nondegenerate with 
level spacing A. The dominant effect of interdot coupling is to introduce a tunneling matrix 
element ta between equivalent single-particle states in nearest-neighbor dots. We neglect 
the tunneling matrix elements between nonequivalent states; this is the usual tight-binding 
approximation, and is justified for nearly identical dots provided is not too large. The 
Hamiltonian is 

H = -Y.(t^ cLcj« + H.C.) +Y.^o. + Y ^ ^j^^j ~ 

(iJ) j,a j 

a 

where i, j are vectors of integers labeling the positions of the dots, cj^ is the creation operator 
for an electron in state a of the jth dot, nj = cj^Cja, and the sum over (i, j) is over nearest 
neighbors only. In the following, we consider arrays of 4 and 6 quantum dots with M = 2, 
3, and 4 single-particle energy levels per dot. 

In this Letter, we calculate the equilibrium electron addition spectrum of (|I|), d{N) /dji = 
kT d"^ In Z/dfi"^, by diagonalizing H and evaluating the grand partition function Z = 
Tr{exp[— (i^' — fiN)/kT]}, where iV = X]j ? A* is the chemical potential, and ( ) denotes 
the thermal average. For the largest array considered, Z involved a sum over 65,536 states. 
Assuming that the tunneling rate from the backgate is much less than the temperature 
kT, the differential self-capacitance of the experimental quantum dot system is given by 
dQ/dV = e^d{N) / dji, which is directly measured by the technique of single-electron capac- 
itance spectroscopy |Tl||T^. 

We first consider an array of four quantum dots arranged in a square with three single- 
particle energy levels per dot, and study the evolution of the addition spectrum as the 
interdot tunneling is increased (see Fig. |I|). Here we have set the tunneling matrix elements 
for all three quantum levels equal to t and have taken A = 0.3?7 and kT = .04U. At 
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t = 0, the behavior characteristic of isolated dots is evident in Fig. ^ namely, three 
peaks separated by the Coulomb blockade energy U + A. Each peak at t = represents the 
addition of four electrons to the array (one to each quantum dot). As t is increased, the 
quantum dot states are split into minibands and each capacitance peak is split into three 
peaks, the central peak representing the addition of two electrons at the same value of the 
chemical potential. When t ~ A, this degeneracy is lifted and one can see three Hubbard 
minibands, each composed of four states. Finally, when t ~ f//2 the energy gap between 
the minibands is no longer discernible. 

Figure |l] shows evidence of three distinct phases in the quantum dot array, which we will 
now analyze in detail. In the weak-tunneling phase, characterized by the degeneracies in 
the addition spectrum, d{N)/dfi has peaks at zero temperature at /i = —2t, 0(x2), and 2t, 
with this pattern repeated, centered at U + A and at 2{U + A). The energies of the lowest 
miniband are those of noninteracting electrons in a tight-binding band with four lattice sites, 
E = —2tcosk, k = 0, ±7r/2,7r. This reflects the fact that at sufficiently small t there is no 
admixture of the higher single-particle states £2 and £3 in the many-body ground state, so 
that the Pauli principle prevents two or more electrons from occupying the same dot, thus 
negating the interaction term in Eq. (Jip. The minibands centered at U + A and 2{U + A) 
are identical, but occur when the states ei and £2 are completely filled. This weak-tunneling 
phase is analogous to the ferromagnetic phase of the Hubbard model in a strong magnetic 
field, the level splitting A playing the role of the field. 

As if: is increased, it becomes energetically favorable to admix the higher single-particle 
states £2 and £3 in the many-body ground state — thereby allowing multiple occupancy of the 
quantum dots — in order to lower the kinetic energy of the system. The degeneracies present 
in the weak-tunneling phase are then lifted by the interaction term in Eq. (|1]). In Fig. |T], this 
transition occurs when t ~ A. For the case of an infinite ID array of quantum dots with M = 
2 levels per dot, this phase transition is equivalent to the ferromagnetic /ant if erromagnetic 
phase transition of the ID Hubbard model in a magnetic field and in the limit U ^ A 
the critical value of the tunneling matrix element is given by |^ 
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where n < 1 is the mean number of electrons per dot. We find that this transition is 
quahtatively similar for M = 3 and 4, and in finite arrays (although it is of course not 
a true phase transition in a finite array). This transition represents a breakdown of the 
Coulomb blockade of individual quantum dots: when t < tc, there is no multiple occupancy 
of the quantum dots in the lowest miniband, while for t > tc, there is a nonzero amplitude 
for multiply occupied dots in the lowest miniband. (It is crucial that the Ea not be spin- 
degenerate, otherwise tc = 0.) Prediction of this observable "phase transition" in quantum 
dot arrays is one of our important new results. 

Despite the destruction of Coulomb blockade in the individual quantum dots when t > tc, 
there is still an energy gap between the minibands above this transition. This energy gap 
is a collective effect which we refer to as "collective Coulomb blockade" (CCB), and is 
analogous to the energy gap in a Mott insulator The CCB regime is characterized 

by strong interdot correlations which are analogous to the antiferromagnetic correlations in 



Mott insulators |]T3[; i.e., occupancy of the state in one quantum dot is anticorrelated 
with the occupancy of that state in its nearest neighbors. 

When t is increased still further, the energy gap between minibands in Fig. |1] collapses. 



The breakdown of CCB ||T^ in the strong-tunneling regime is directly analogous to the 
Mott-Hubbard insulator-metal transition For M = 2 single-particle energy levels per 

dot, this transition occurs at t/U = oo for a ID array |T3[, but for M > 2, the transition 
is expected to occur at a finite value of t due to the absence of Fermi-surface nesting [0,|l6 



Fig. ID is consistent with the critical value obtained in Ref. for an infinite array with 



M = 3, tc/U = 0.39 {although the metal-insulator transition studied in Ref. ||T6 



occurs m 



an unspecified model which is merely analogous to Eq. (^}. The signature of this "metal- 
insulator" transition in the electron addition spectrum is relatively insensitive to finite-size 
effects (which can have a drastic effect on the conductivity of the system ||^). 

Figure ^ shows the effects of disorder and magnetic field on the degeneracies which occur 



in the addition spectrum in the weak-tunnehng regime. d{N) / dfi is plotted in gray scale, 
where white represents the peaks. The data are for the lowest miniband of a singly-connected 
array of six quantum dots with two single-particle energy levels per dot. The area enclosed 
by the "ring" of dots is (llOnm)^, and the system parameters are U = lOmV, A = 3mV, 
t = 1.5mV, and T = I.IK. Disorder with a magnitude ~ 0.2mV is included in the on-site 
energies, and one of the t's is taken to have a negative sign in zero field. The magnetic field 
modifies the interdot tunneling matrix elements in Eq. (|1|) by a Peierls phase factor (we 
neglect modifications to the magnitude of t). The modification of the quantum confinement 
by the magnetic field may be taken into account [0 by setting ei = [A^ + {ujc/'^y]^^'^ — A, 
€2 = 2[A^ + — A — uJc/2, where Uc is the cyclotron energy (the effective mass 

of bulk GaAs has been used). The Zeeman term is only .09mV at lOT, and is neglected. 
Disorder mixes the states in the miniband with ±k and splits the degeneracies present at 
zero field by an amount ~ (ef — e^^)^/^ (disorder in the hopping matrix elements has a similar 
effect). This splitting is obscured in Fig. 0by the finite electron temperature. The magnetic 
field has two quite different effects: it explicitly breaks time-reversal invariance, leading to 
a splitting proportional to the field; in addition, because the splitting between ei and £2 is a 
decreasing function of magnetic field, the magnetic field can actually drive the system into 
the CCB regime, leading to an abrupt splitting. In Fig. ^, this field-induced phase transition 
is clearly visible in the highest doublet at ~ IT, and in the intermediate doublet at ~ 2T, 
but is absent in the lowest doublet, in agreement with the density dependence of Eq. (^. 

The experimental capacitance spectrum obtained by Ashoori et al. [|TT1 for a disordered 
quantum dot with several local minima — which can be thought of as a random array of 
small quantum dots between 130A and 360A in size — exhibits degeneracies much like those 
in Fig. ^ which are split by a magnetic field. If these degeneracies are to be ascribed 
to the above mechanism, it would require that the local minima be arranged in a nearly 
symmetrical pattern (it is not necessary that they form a ring — similar degeneracies occur in 
other symmetrical arrays.) This possibility is consistent with the fact that the degeneracies 
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disappeared after the sample was thermally cycled [|TT|. Of course, it is crucial for this 



interpretation that the energy levels in the local minima be non-spin-degenerate, which 
would require spin-orbit coupling, or some other mechanism. To our knowledge, no other 
explanation has been advanced for these anamolous degeneracies. 

Finally, we consider the effects of energy-dependent tunneling. To this point, we have 
assumed equal tunneling matrix elements for each of the quantum dot states, but in fact 
the more weakly bound states should be connected by stronger tunneling matrix elements. 
To simulate this effect, in Fig. § we show the electron addition spectrum for a linear array 
of four quantum dots with four single-particle energy levels per dot with tunneling matrix 
elements ta = 0.1f/(1.5)°~^, a = 1, ... ,4. The level spacing is A = 0.3U. Note that for a 
linear array with open boundary conditions, the degeneracies discussed above do not occur. 
When disorder was introduced into the tunneling matrix elements, the gross features of 
the spectrum were unchanged, but the peak spacings within the minibands were different. 
The energy gap between the minibands decreases with increasing chemical potential; when 
the value of U was decreased 40% relative to A and ta, the energy gap between the third 
and fourth minibands was no longer discernible. As discussed above, this breakdown of 



CCB is analogous to the Mott-Hubbard insulator-metal transition Quantum dot 

arrays should thus provide a system where this transition can be studied by tuning a gate 
voltage, which would be inconceivable in the metal-oxide systems which display conventional 
Mott-Hubbard metal-insulator transitions. 

A reexamination of the work of Kouwenhoven et al. [0] leads us to conclude that the 
analogue of the Mott-Hubbard metal-insulator transition may have already been observed 
in a ID array of quantum dots. In Ref. 0, the conductance spectrum of a linear array 
of fifteen GaAs quantum dots ~ lOOnm in diameter was measured. While our capacitance 
spectra cannot be compared directly to those conductance spectra, the peak positions should 
be the same in linear response, though the peak heights and possibly the line shapes would be 
different. At low values of the gate voltage Vg2 (corresponding to low values of the chemical 
potential in the array), conductance peaks separated by 4mV were seen in Ref. |0. Given 



the electrostatic "lever arm" of the gate with respect to the quantum dots of about 4, this 
implies a ImV spacing in energy between the peaks, which is roughly the Coulomb blockade 
energy for a lOOnm GaAs quantum dot. We therefore interpret these peaks as Hubbard 
minibands in which the individual discrete states are unresolved due to a finite electron 
temperature (c./., the lowest miniband in Fig. As the gate voltage Vg2 was increased, 
the conductance peaks were seen to broaden 0, and some structure in the peaks became 
observable; finally the gap structure with a period of 4mV was destroyed, which we interpret 
as the breakdown of CCB. The similarity to Fig. § is striking. Further experimental work 
is clearly warranted to map out more fully the "phase diagram" of quantum dot arrays. 

Some shortcomings of the model considered here [Eq. (|l])] should be pointed out. The 
tight-binding approximation employed here to treat interdot coupling is probably not ad- 
equate to treat the strong-tunneling regime where the Coulomb blockade energy gap has 
broken down completely (the plateau region in the conductance spectra of Ref. 0); the 
analysis given in Refs. |^,^ of that regime is probably more appropriate. The Hubbard- 
type interaction term in Eq. (|I]) is also a rather severe idealization of the electron-electron 
interactions within the quantum dots. More realistic treatments of the electron-electron in- 
teractions within a quantum dot are necessary [|T^] to treat the regime of strong magnetic 



fields, where strongly-correlated precursors of the bulk-incompressible fractional quantum 
Hall states are expected to occur. However, in the regime of weak tunneling and weak mag- 
netic fields which we are primarily concerned with, we believe the collective phenomena in 
semiconductor quantum dot arrays are adequately described by our model. 

In conclusion, we have calculated the electron addition spectrum, or differential capac- 
itance, of arrays of 4 to 6 coupled semiconductor quantum dots in and near the Coulomb 
blockade regime using a generalized Hubbard model. We argue on general grounds that the 
interplay of quantum confinement, interdot tunneling, and strong intradot Coulomb interac- 
tions leads to three distinct zero temperature phases in quantum dot arrays which have clear 
signatures in the electron addition spectrum: for weak tunneling, the quantum dot states are 
split into minibands but the Coulomb blockade of individual dots is maintained; for inter- 
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mediate tunneling strengths, the Coulomb blockade of individual dots is destroyed but there 
remains an energy gap of collective origin analogous to the energy gap in a Mott insulator; 
and for strong tunneling, the Coulomb blockade is destroyed altogether. The weak-tunneling 
phase is characterized by degeneracies in the electron addition spectrum which are lifted by 
disorder and by magnetic fields. The intermediate phase we refer to as "collective Coulomb 
blockade." The breakdown of collective Coulomb blockade in the strong tunneling regime is 
analogous to the Mott-Hubbard insulator- metal transition We have argued that the 

phenomenology of an early experiment on transport through a ID array of quantum dots 
0] can be understood in this framework. We believe that quantum dot arrays are rather 
ideal physical systems to directly verify the predictions of Hubbard-type strongly-correlated 
electron interaction models, albeit at very different energy scales (~ meV) compared with 
what these models were originally intended for (~ eV). 

This work was supported by the National Science Foundation (NSF) and the United 
States Office of Naval Research (US-ONR). 
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FIGURES 

FIG. 1. A plot showing the evolution of the electron addition spectrum d{N)/dn, as a function of 
the interdot tunneling matrix element t for an array of four quantum dots with three single-particle 
energy levels per dot, arranged in a square. 

FIG. 2. A gray scale image of the electron addition spectrum versus magnetic field in the 
lowest miniband of a ring of six quantum dots with two single-particle energy levels per dot. Here 
U = lOmV, A = 3mV, t = 1.5mV, and T = I.IK. The radius of the ring is 62nm and the effective 
mass for bulk GaAs has been used. 

FIG. 3. A plot of the electron addition spectrum d{N)/diJ, versus the chemical potential fi for a 
linear array of four quantum dots with four single-particle energy levels per dot. The level spacing 
is A = 0.3f7, and the tunneling matrix elements are ta = 0.1J7(1.5)"~^, a = 1, ... ,4. Solid curve: 
kT = mU. Dotted curve: kT = .02U. 



12 





2 4 6 8 10 
Magnetic Field (T) 



